Whispering gallery mode enhanced optical force with resonant tunneling excitation in the 

Kretschmann geometry 
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The boundary element method is applied to investigate the optical forces when whispering gallery modes 
(WGMs) are excited by a total internally reflected wave. Such evanescent wave is particularly effective in 
exciting the high-Q WGM, while the low angular or high radial order modes are suppressed relatively. This 
results in a large contrast between the forces on and off resonance, and thus allows for high size-selectivity. We 
fully incorporate the prism-particle interaction and found that the optical force behaves differently at different 
separations. Optimal separation is found which corresponds to a compromise between intensity and Q factor. 



Microcavities with whispering gallery mode (WGM) are 
important components in nanophotonics, nonlinear optics, 
quantum optics, and many other areas. 1 As it is very dif- 
ficult to accurately control the resonance frequency or the 
size of a micro-cavity during its fabrication, a convenient 
way to pick up the particles with the desired resonance fre- 
quencies or size is highly desirable. Recently, it was pro- 
posed in Ref. 2 that the optical force acting on a microsphere 
under an evanescent wave illumination 3,4,5 ' 6 ' 7 is highly size- 
selective: at resonance the microsphere experiences a siz- 
able force, whereas at off resonance the force is negligible. 8 
Utilizing this property and the high quality factor (Q > 10 7 ) 
of WGM, 1,9 the size-selective force would allow for accu- 
rate size-sorting of microsphere cavities, with an accuracy 
of ~ 1/Q. However, the generation of evanescent wave re- 
quires a substrate. Ref. 2 considers only the regime where 
the microsphere is away from the near field of the substrate, 
and the interaction of the microsphere and substrate is only 
treated in a phenomenological manner. 

In this article, using the boundary element method 
(BEM) 10 and the Maxwell stress tensor, 11 we intend to go 
beyond Ref. 2 by fully incorporating the influence of the 
substrate. We confine ourselves to 2D calculations with 
a transverse electric (TE) polarization, since 3D calcula- 
tion is difficult due to the need of extremely fine meshes 
in the BEM calculation. We expect that our result should 
be qualitatively relevant for the 3D system. We consider the 
Kretschmann geometry as depicted in Fig. 1(f): a micro- 
cylinder (MC) of radius a sits on top of a prism, both have a 
refractive index of n p = 1.5. The combined system (in air) 
is then illuminated by an incident plane wave propagating at 
6 = 45°. As the critical angle of the prism is 6 C = 41.8°, 
the incident wave is total internally reflected at the prism 
surface, and an evanescent wave is generated there. The 
transfer of photon momentum to the MC is allowed by the 
scattering or tunneling of the evanescent wave. We shall see 
that when the particle is not in the near field of the prism 
surface, the treatment of Ref. 2 gives a correct description 
of the situation. For the near field case, we observe a series 
of interesting phenomena, some are caused by the interplay 



between the increasing intensity and the decreasing Q as the 
particle approaches the surface, and some are caused by the 
interaction between the particle and the substrate. 

The research on WGM-induced forces 12 ' 13,14,15 ' 16,17,18 can 
be dated back to the early eighties when Ashkin reported 
the first observation of the force in air. 19 Later, the WGM- 
induced forces were also observed by Fontes et al in wa- 
ter. 20 Nevertheless, only a modest enhancement in the opti- 
cal force has been achieved owing to the use of propagat- 
ing waves in these works. Meanwhile, it is well-known 
that the Kretschmann geometry can efficiently excite the 
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FIG. 1 : (Color online) Optical forces at different height from the 
prism as functions of the wavelength. (Blue) dotted line: a micro- 
cylinder illuminated by a plane wave in the absence of the prism. 
(Black) solid line: the horizontal force F x acting on a micro- 
cylinder with the prism. (Red) dashed line: the vertical force F y 
with the prism, (a) d — 1.1a, (b) d = 1.2a, (c) d = 1.3a, Inset 
of (c): an enlarged section of (c). (d) d = 1.1a and (e) d = 1.2a 
for shorter wavelength. The Qj so and resonant wavelength for the 
typical modes in (d) and (e) are, respectively, 5399 and 0.31110a 
for TE26, 7632 and 0.30047a for TE^ 7 , 10817 and 0.29055a for 
TE^g, 15368 and 0.28128a for TE^ 9 , 96 and 0.31044a for TEl 2 , 
115 and 0.29928a for TEl 3 , 139 and 0.28892a for TEl 4 , and 170 
and 0.27926a for TE25. (f) Illustration of the Kretschmann ge- 
ometry, where the evanescent wave is generated by total internal 
reflection on top of the prism. 
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WGM of a high-Q microsphere. Kawata and Sugiura 3 
demonstrated versatile optical manipulation by evanescent 
waves subsequently. 21 Although the Kretschmann geometry 
is well-known to be suitable for both optical manipulation 
and the excitation of WGM, the resonant forces associated 
with evanescent wave excited WGM have never been con- 
sidered experimentally. Here, we hope to provide a more 
complete treatment for the WGM induced forces, and hope 
our work can stimulate experimental works in this area. 

In Figs. l(a)-l(c), we show the optical forces acting on a 
MC for A = 0.32a - 0.4a (size parameter k p a « 24 - 30). 
The (blue) dotted line represents both the horizontal and ver- 
tical forces (equal by symmetry) acting on a bare MC with- 
out the prism. 22 The (black) solid line represents the hori- 
zontal force F x and the (red) dashed line represents the ver- 
tical force F y , in the presence of the prism. For the case 
without the prism, the plane incident wave hits directly on 
the MC, and for the case with the prism, the plane wave is 
partially reflected by the prism surfaces and partially con- 
verted into an evanescent wave on the prism's top surface. 
In the case of a bare MC, it can be seen from Fig. 1(a) that 
the force at resonance is comparable to the force at off res- 
onance; thus size-selectivity is weak. On the other hand, 
the case with an incident evanescent wave (generated by the 
prism) is very size-selective, as we shall see below. 

Let us first consider the horizontal force F x . As shown 
in Fig. 1(c) where the separation d = 1.3a, it is quite clear 
that the force on first radial order resonance is significantly 
stronger than the force at off resonance, and force on the 
second order resonances can hardly be seen at all. This re- 
sults in a high peak-to-baseline ratio, defined as the peak 
resonance force divided by the maximum of off resonance 
or lower-Q resonance force. This would imply the possibil- 
ity of size-selective manipulation: upon illumination, only 
those particles whose sizes happen to be in resonance with 
the incident wave will experience a sizable optical force, 
whereas the off resonance particle will be left untouched. In 
the limit of large separation, the peak-to-baseline ratio can 
be several orders of magnitude higher than what is predicted 
here, as found out in Ref. 2. Indeed, the results presented in 
Fig. 1(c) already begin to resemble that of Ref. 2. 

Although it is highly desirable to have a high peak-to- 
baseline ratio, it is also important to achieve large force. 
This case is illustrated in Figs. 1(a) and 1(b). As the sep- 
aration decreases from d = 1.3a in Fig. 1(c) to d = 1.2a 
in Fig. 1(b), the force on first order resonance increases by 
a factor of ~2-3 as the intensity becomes higher in the near 
field, while its linewidth is broadened slightly. At d = 1.2a, 
the second order resonance can now be observed, due to 
the higher intensity. Although the peak-to-baseline ratio for 
d = 1.2a is smaller than that of d = 1.3a, the size selectivity 
is still very strong. At an even closer separation [d = 1.1a 
in Fig. 1(a)], although the intensity of the evanescent wave 
increases, the strength of the first order resonance force re- 
mains at the same level as when d = 1.2a. This is a conse- 
quence of the Lorentzian lineshape, 2 which implies that the 
resonance force is inversely proportional to the linewidth, or 
directly proportional to its Q. As d decreases, the stronger 
modal coupling between the prism and the MC broadens 
the linewidth (or decrease in Q), and causes the resonance 
force to decrease, compensating the increase in the local in- 
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FIG. 2: (Color online) Optical forces as a function of cylinder- 
prism separation d at various situations, (a) at the TE22 resonance, 
(b) at an off resonance frequency between TE22 and TE^g, (c) at 
the TEi 8 resonance. 

tensity I\ oc . The background force and the forces associated 
with the second order resonance have increased consider- 
ably, making the peak-to-baseline ratio small. We note that 
as the MC approaches the prism, the resonance frequencies 
do not shift and it is the Q's that drops. This decrease in Q 
will significantly reduce the size-sensitivity of the resonance 
force, which scales like ~1/Q. 

A more explicit demonstration on the selectivity due to 
high-Q is shown in Fig. 1(d) and 1(e) which show the op- 
tical forces with the presence of the prism in higher size 
parameter regime A = 0.28a — 0.32a (fc p a w 30 — 34). It is 
seen in Fig. 1(d) that with d = 1.1a several radially second 
order WGMs can be well excited, resulting in comparable 
optical forces of the radially first order WGMs as indicated 
by the downward arrows. We find that the Q value of the 
TE23 mode exceeds 10 2 (see the details in the caption of 
Fig. 1). The MC with such a value of Q can already collect 
sufficient incident field energy that yields a non-negligible 
amount of light scattering. Figure 1(e) shows that for in- 
creased separation to d = 1.2a, these second order modes 
remain not well excited again, simply because of the de- 
creased local field. At this separation, however, the radially 
first order WGMs with Q in the order of 10 3 and 10 4 can be 
well excited and are manifested as the several sharp peaks 
[as labeled in Fig. 1(e)]. 

Figure 2 clarifies the separation dependence of the opti- 
cal force on and off resonance. Here we have chosen three 
different frequencies, one for TE22 (with quality factor of 
isolated particle <3j so = 1395) shown in Fig. 2(a), one for 
TE 2 8 (Q iso = 49) 10 shown in Fig. 2(c), and the last one for 
an off resonance frequency shown in Fig. 2(b). We remark 
that the strength of the resonance force depends on both 
the local intensity 7i oc and the mode's Q, i.e. F oc QI\ oc - 
The observed Q in the force spectrum is in fact related to 
the original quality factor for the isolated MC, Qj so , and 
the coupling Q-factor, Qcoup, through the equationl/Q = 
1/Qiso + l/<2coup- ln Fig. 2(b) and 2(c), the optical force 
decreases monotonically as d increases. The off resonance 
case of Fig. 2(b) is straightforward to understand: as d in- 
creases, the local intensity decreases, so the force decreases 
as well. For the second order resonanceTE 2 g case plotted 
in Fig. 2(c), the system's Q is completely determined by its 
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FIG. 3: (Color online) Near fields for high-Q (left column) and 
low-Q (right column) resonances, (a) and (c) are, respectively, 
for TE22 and TEf g in the absence of the prism, (b) and (d) are, 
respectively, for the same resonances as (a) and (c), but with the 
prism at a separation d — 1.2a (lower panels). 



small Qi so = 49 (because Q; so -c Qcoup), consequently as 
d increases, the total Q is essentially unchanged while ii oc 
decreases, causing the optical force to drop. For the high-Q 
resonance of TE^ shown in Fig. 2(a), it can be anticipated 
that the optical forces will be weak when the MC is either 
too close to (due to a small Q) or far away from (due to 
small ii oc ) the prism, and there is a peaked force in between. 
In Fig. 2(a), this peak occurs at d w 1.16a. 

To give an overall perspective of the effect of Q value and 
the effect of the prism, we compare the near fields in Fig. 3 
for four different situations corresponding to Figs. 2(a) and 
2(c). These figures clearly demonstrate the resonant tunnel- 
ing mechanism which accounts for the optical forces. We 
can draw a conclusion that a high Q value is the prerequisite 
for the resonant tunneling. Therefore in the scheme with the 
prism for optical sorting operated at a specific frequency, 
high Q particles can be easily distinguished and selected. 

We finally consider the vertical force F y . As Rayleigh 
dielectric objects tend to move to region of high intensity 
in order to minimize the electromagnetic free energy, one 
may believe that the vertical force is attractive relative to 
the prism. But as shown in Fig. 1(c) where d = 1.3a, the 
vertical force versus wavelength curve actually consists of 
a background repulsion which lifts the entire curve upward 
plus a peak and dip pair at each resonance position. As the 
MC approaches the prism [Figs. 1(a) and 1(b)], the line- 



shape of F y evolves and deviates from that of the original 
lineshape in Fig. 1(c). We see that the resonance force has 
become stronger, and F y can now be either attractive (at the 
dip) or repulsive (at the peak). This resonance lineshape 
(the peak and dip pair) can be understood as follows: when 
the frequency changes across the resonance frequency, the 
WGM switches from in phase to out of phase with the in- 
cident field, and this results in a change of the sign in the 
force. We note that in addition to typical resonance peak- 
and-dip, there is still a bias of the force to the repulsive side. 
Heuristically, this repulsive bias can be traced to the fact that 
the incident photon momentum is upward, so it should not 
be surprising to see the photons pushing the MC up. We 
note that in the 1992 experiment of Kawata and Sugiura, 3 a 
repulsive force is also observed. This subtle effect warrants 
further studies. 

In summary, we have numerically investigated the WGM 
enhanced optical force as excited by an evanescent wave. 
The high peak-to-baseline ratio associated with evanescent 
wave excitation makes it highly size-selective and may thus 
be useful for particle sorting. The notion of Ref. 2 is con- 
firmed even in the presence of a substrate, and we have gone 
beyond Ref. 2 and show that the idea even works in the sub- 
strate's near field. In general, a high Q mode guarantees 
a stronger optical force, as it favors the resonant tunneling 
of the evanescent light. As a consequence of the competi- 
tion between the increased coupling loss at small separation 
and the decreased intensity at large separation, there exists 
an optimal particle-prism separation in which the resonance 
force is maximized. Accordingly, in an optical sorting ex- 
periment, it would be highly desirable to have some other 
means to keep the particles at an appropriate height, for 
if the MC is too close to the prism surface, the linewidth 
broadening would reduce the sorting accuracy as well as the 
peak-to-baseline ratio, but if the MC is too far away, the 
force is too weak. Finally, we remark that the induced op- 
tical force can further be enhanced by coating the substrate 
with a metallic coating that supports surface plasmon, 23 or 
by coating the prism with a dielectric cavity layer. 24 
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